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Abstract — We report on charge transport and current fluc- 
tuations in a single bacteriorhodpsin protein in a wide range 
of applied voltages covering direct and injection tunnelling 
regimes. The satisfactory agreement between theory and available 
experiments validates the physical plausibility of the model 
developed here. In particular, we predict a rather abrupt increase 
of the variance of current fluctuations in concomitance with that 
of the I-V characteristic. The sharp increase, for about five orders 
of magnitude of current variance is associated with the opening of 
low resistance paths responsible for the sharp increase of the I-V 
characteristics. A strong non-Gaussian behavior of the associated 
probability distribution function is further detected by numerical 
calculations. 

I. Introduction 

One of the more intriguing properties of proteins is the 
strong interdependence of function and structure. Each pro- 
tein consists of a specific sequence of aminoacids, spatially 
organized (folded) in few possible structures, each of them 
corresponding to a specific level of functioning. In the un- 
folded state the protein does not work. The switching from 
a structure to an other (say, for example, from the native to 
the active state) allows the protein to change its activity. This 
feature has numerous technological applications, therefore it 
has recently received much attention from different fields of 
investigation. In particular, the structure switching can be used 
in the passive way to reveal the changing of functionality as 
induced by external agents. This is the case, for example, of 
recent investigations that point to the conformational change 
for revealing the capture of an external agent, like an odour 

m. 

The structure switching could also be used in the active 
way, i.e. by inducing it with an external agent to modify some 
peculiar protein properties. An example of this application 
was obtained by a group at the Weizmann Institute [2| which 
using a light sensitive protein, the bacteriorhodopsin (bR), 
was able to reveal a significant increase of the electrical 
current when the protein is irradiated with green light. This 
particular feature exploits the possibility to use the bR as the 
active element of an organic field effect transistor (OFET) 
0, with the interesting traits of low cost, easy production 
and efficiency on a large range of voltages. In particular, the 
protein exhibits a good stability under the application of high 
fields, showing a response which explores the range from 
direct to injection (Fowler-Nordheim) tunneling regime. The 
interpretation and the theoretical modeling of this response 



was previously carried out by neglecting the protein internal 
structure (4) and therefore it was not able to describe the 
modifications induced by the conformational change. 

In this paper we propose a model of the electrical properties 
of proteins, based on the assumption that their modifications 
follow by their conformational change. The model has a wide 
applicability and allows the exploration of the topological 
properties of the proteins, such as of the static and dynamical 
responses under an electrical solicitations |5|. Therefore, we 
are able both to describe the existing data and to provide some 
expectations. In particular, we reproduce the I-V characteristic 
given by experiments [6| and also predict the associated 
fluctuations of the steady current. Finally, by analyzing the 
distributions of current fluctuations for different voltages, we 
find that they follow the universal behaviour originally found 
by Bramwell-Holdsworth-Pintor (BHP) [7|, and continuously 
rediscovered in most critical system [7], [8|, @ 

II. Experiment and model 

Bacteriorhodopsin is contained in a part of the cellular 
membrane, the purple membrane (PM), of the halophile 
Halobacterium salinarum ifTOl . PM is a particularly simple 
kind of cellular membrane, being constituted only by trimers 
of bR and some lipids to stabilize them, and is quite easy to 
produce. It appears as a quasi-2D structure, with a thickness 
of about 5 nm, corresponding to the protein dimension. These 
characteristics make it very suitable for any kind of experi- 
mental investigations. 

A seminal set of papers [2| reported the measurements car- 
ried on metal-insulator-metal (MIM) junctions of millimetric 
diameter, with the insulator constituted by a 5 nm monolayer 
of PM. The current response (nA level) was investigated in a 
small range of bias (0-^1 V), and in this range it appeared as 
super-Ohmic, growing overall a factor of 2 when the sample 
is irradiated by a green light. These results suggest that in this 
protein, like in some organic polymers [1 1J, electrical transport 
(ET) is ruled by tunneling mechanisms and that it is strongly 
related to the protein structure fl2l . In successive experiments 
fl3l . the I-V characterization was performed at the nanometric 
scale, with the conductive atomic force microscopy (c-AFM) 
technique. In this technique, one of the contacts is constituted 
by the tip of the c-AFM (100 -f- 200 nm of nominal radius). 
With respect to previous measurements, in the common bias 
range, the current response is lower for about 4 order of 



magnitude, thus the sample is able to sustain higher voltage 
up to about 5 4- 10 V. Measurements were performed without 
any extra light irradiation. At voltages above about 2.5 V, 
the presence of a cross-over between the direct tunneling 
regime and the injection or Fowler-Nordheim (FN) tunneling 
regime lTT4l was evidenced. To our knowledge, a microscopic 
interpretation of the FN tunneling accounting for the structure 
of the protein is not present in literature, while some macro- 
scopic approaches, which deal the protein as a uniform bulk 
insulator can be found |4|. This paper aims to fill this gap 
by giving a tool for modeling the proteins at a microscopic 
level, and, at the same time, able to reproduce the macroscopic 
data. The model, hereafter called INPA (impedance network 
protein analogue), is based on a coarse-grained description 
of the protein: The aminoacids, taken as single centers of 
interaction, are the unique responsible of charge transfer and/or 
polarization. The protein is represented with an impedance 
network, which retains the main topological features of the 
protein lfT2l . with a mechanism of electrical response that 
depends on the kind of applied bias (DC, AC) 0. In the 
present case, by using the analogy with the charge transport 
among localized states, we model the charge transport through 
a single protein by means of a sequential tunneling mechanism 

The approach integrates structure and function of the protein 
and, therefore, it constitutes a fundamental step toward the 
use of proteins like PM as an active part in bio-devices. With 
respect to previous approaches, INPA has several advantages: 
i) it can explain the different electrical responses at increasing 
voltages and in the presence/absence of green light; u) it can 
reconcile the interpretation of different data 12), fl3l : Hi) it 
can be applied to other proteins whose tertiary structure is 
known. 

In brief the layout of the INPA is as follows. The protein is 
mapped into an impedance network, by using the C Q atoms of 
each aminoacid as nodes. The number of links for each node 
is determined by the cut-off interaction radius, say Rq- In the 
present analysis, we choose Rq = 6 A, a value that optimizes 
the native to activated state resolution 1121 . Each link is 
associated with an impedance (a simple resistance in this case) 
whose value depends on the distance between amino-acids as: 
r i,j = Pkj/Aij, where p indicates the resistivity, here taken 
to be same for all the links, which in general depends on the 
voltage as detailed below, the pedices i,j refer to the amino- 
acids between which the link is stretched, lij is the distance 
between the labeled aminoacids taken as point like centers and 
Ai.j is the cross-sectional area shared by the labeled amino- 
acids: Ai d =ir(E&-t$j/4). 

To take into account the superlinear I-V response, the link 
resistivity, p, is chosen to depend on the voltage drop of each 
resistance as: 

( Pmax eV < $ 

P(V) = { 

{ pMAx{£r) + p mm (l - eV > $ 

(1) 

where pmax is the resistivity value which should be used to 



fit the I-V characteristic at the lowest voltages, p m i n <C Pmax 
plays the role of an extremely low series resistance, limiting 
the current at the highest voltages, and $ is the threshold en- 
ergy separating the two tunneling regimes (a kind of effective 
height of a tunneling barrier). Since ET is here interpreted in 
terms of a sequential tunneling between neighboring amino- 
acids, the above interpolation formula reflects the different 
voltage dependence in the prefactor of the current expression 
(4): / ~ V in the direct tunneling regime, and I ~ V 2 in the 
FN tunneling regime. 

For the transmission probability of the tunneling mechanism 
we take the expression given by Ref. Ifl4ll : 
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(3) 

where Vj j is the local potential drop between the couple of 
i,j amino-acids and m is the electron effective mass, here 
taken the same of the bare value. To model bacteriorhodopsin 
in its native state (in dark) lfT2l . lfT31l we have taken the PDB 
entry 2NTU. The network is then studied with a point-like 
contacts configuration, connected to the external bias by means 
of perfectly conductive contacts, the input on the first amino- 
acid and the output on the last amino-acid of the primary 
structure. 

III. Results 

Figure ([T]) reports the current evolution in bR as a function 
of time (measured in number of iteration steps) for three 
different applied voltages: 1 V, in the low linear region; V=3.5 
V in the transition region; V=5 V, approaching the extreme 
linear regime. The data relative to 5 V have been shifted up 
of 100 nA in order to better resolve the figure. The current 
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Fig. 1. Current signal vs time (measured in number of iteration steps) at 
three different voltage values, from bottom to top: 1 V, 3.5 V, 5 V. 

response exhibits a low level of noise at low voltages, where 



direct tunneling is the dominant transport process. Near the 
transition between the two tunneling regimes, the amplitude 
of fluctuations increases, and finally, at high voltages an 
increasing number of very high current events (up to five 
order of magnitude larger than the major part of the signal) is 
manifested. 

Figure Q compares the numerical and experimental data 
of the current voltage characteristic obtained with an AFM 
technique carried out on a bacteriorhodopsin monolayer at 
room temperature lfl3ll . As relevant computational inputs, 
we have taken: $, the threshold voltage value, of 219 meV, 
p MAX = 4 • 10 13 Q A, p mm = 4 • 10 5 f2 A. The satisfac- 
tory agreement between theory and experiments validates the 
physical plausibility of the model developed here. The sharp 
increase of the current above about 3V is associated with the 
cross-over between direct and injection sequential-tunneling 
regimes. 

Figure ^ reports the variance of current fluctuations cor- 
responding to the I-V characteristic of Fig. Here we 
notice a rather abrupt increase of the variance of current 
fluctuations in concomitance with the cross-over region of the 
I-V characteristic. The sharp increase for about five orders of 
magnitude of current variance is associated with the opening of 
low resistance paths responsible for the sharp increase of the 
I-V characteristics. By analising the frequency distribution 
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Fig. 2. I-V characteristic for bR at the nanometric scale: Experiments (bold 
line) 1131 vs. simulations (circles), tiny line is a guide to eyes. 

of currents, v, for different voltages, a strong non-Gaussian 
behaviour can be observed (see Fig.Q). This result is not 
completely surprising, being partially due to the small size 
of the system under consideration |9| but also describing the 
critical condition of a system in which, for all the considered 
voltages, two different regimes are in competition. On the 
other hand, this kind of systems, and their fluctuations in 
particular, have received in the last few years Q a huge 
attention, mainly due to the possibility to reconduct them 
to a 'universal' behaviour. This behaviour is described by a 
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Fig. 3. Variance of current fluctuations obtained from simulations 
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Fig. 4. Frequency distribution of currents for different values of applied 
voltages. 

generalized Gumbel distribution, one of the reference distribu- 
tions for the statistical modeling of extreme events [16|. This 
distribution was extensively studied by Bramwell-Holdsworfh- 
Pintor (BHP) who found it in very different contexts going 
from fiuidynamics, to self-organized critical systems, and to 
resistance fluctuations. The BHP distribution function has the 
following expression 

p(z)=A{e- z ~^) a (4) 

where p(z) is the probability density function (PDF). We found 
that the skewed distributions reported in Fig.Q can be simply 
recolapsed to a unique curve, by using as variable the natural 
logarithm of current and adopting the rescaled quantities: 

n(z) = a$(x), z = (x- < x >)/a (5) 

with &(x), < x >, a assuming the standard meaning of the 
probability density function, the mean value of the variable 



and the standard deviation, respectively. The choice of the 
logarithm instead of the current is suggested by the slow 
dependence of this quantity on the applied voltage Q, 10. In 
particular, in a large voltage range, the distributions collapse 
and can be described by the expression: 

p(z) = ae -( a * +7) - e ~ (ax+7) (6) 

with a = 1.3 and 7 the Euler constant |8|. 
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Fig. 5. Rescaled probability density function of current fluctuations for 
different values of applied voltage. The bold line is the scaled Gumbel 
distribution with a=1.3 



IV. Conclusion 

The paper reports a microscopic model of current and 
fluctuations in proteins. In particular PM patches of the light 
receptor bacteriorhodopsin have been analyzed in a very wide 
range of voltage (0.0 IV - 7 V) where the current response 
is found to exhibit two different regimes of charge transport, 
direct and FN tunneling. A unique microscopic description, 
the so-called INPA model, is able to take into account both 
these regimes, thus producing the current response with a good 
agreement with available experimental data. Furthermore, the 
current fluctuations are found to follow a universal behaviour 
typical of systems close to a phase transition. 
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